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ABSTRACT 

Gravitational lensing is an important tool for probing the mass distribu- 
tion of galaxies. In this letter we report the discovery of a new quadruple 
lens HST 1411+5211 found in archived WFPC2 images of the galaxy cluster 
CL140933+5226. If the galaxy is a cluster member then its redshift is 2; = 0.46. 
The images of the source appear unresolved in the WFC implying that the source 
is a quasar. We have modeled the lens as both a single galaxy and a galaxy plus 
a cluster. The latter model yields excellent fits to the image positions along with 
reasonable parameters for the galaxy and cluster making HST 1411+5211 a likely 
gravitational lens. Determination of the source redshift and confirmation of the 
lens redshift would allow us to put strong constraints on the mass distribution 
of the lensing galaxy. 

Subject headings: galaxies: fundamental parameters-gravitational lens- 
ing 
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1. Introduction 

Gravitationally lensed systems have a variety of astronomical uses. The most 
obvious is that they provide a direct probe of foreground mass distributions 
(stars, galaxies, and clusters). Systems which exhibit multiple images of a single 
source have an additional application; if the source varies, the system can be 
used to measure Bubble's constant ( |Kefsdal 1964| ). Additionally, the statistics 



of multiple gravitational lenses can provide a direct and powerful test of cosmic 
structure formation theories ([Turner et al 198^) . 



To date eight quadruple lenses have been found, seven with observed lensing 
galaxies. Detailed modeling of these systems, has revealed the need for substan- 
tial external shear in two of the lenses and misaligment between the major axes 
of the light and matter in one of the lenses ( [Witt fc Mao 1997 ). Increasing the 



database of known quadruple lens will improve our knowledge of galaxy mass 
distributions, allowing us to put tight constraints on the mass-to-light ratios and 
hence the dark matter content in the inner few kiloparsecs. If sufficent lenses 
can be found over a large redshift range evolution of galaxy mass-to-light ratios 
could also be measured. 

In this letter we report the discovery of a new quadruple lens system found 
in archived WFPC2 images. The primary lens, which is clearly identified, is a 
galaxy located near and quite likely in a cluster at z = 0.46. In §0 we discuss 
the observations and the optical appearance of the lensing galaxy. In §^ we give 
a brief review of gravitational lensing physics and apply one single-lens and two 
double-lens models to the data. Throughout the paper, we use f2o = 1 and use 
h =Ho/100. 



2. Observations 

Observations of the cluster CL140933+5226 (3c295) were obtained on 25 July 
1994 with the Wide-Field Planetary Camera on HST (images were obtained 
from the archive). The total integration time was 12600 seconds in 6 individual 
exposures using the F702W filter. The 6 images (archive IDs are u2c40a01t, 
u2c40a03t, u2c40a05t, u2c40a02t, u2c40a04t, u2c40a06t.c0h) were registered and 
combined using an average with cosmic ray (positive) outlier rejection. 

The images were reduced and analysed as part of a larger program to per- 
form two-dimensional surface photometry to derive structural and photometric 
parameters for elliptical galaxies in clusters over a wide range of redshift (|Schade| 



et al. 19971) . Each fitted galaxy was inspected for goodness of fit and the presence 



of irregular residuals. The lensed background images were not recognized until 
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the best-fit galaxy model was subtracted (see Figure 1). 

The lensing galaxy is Number 162 from Table 6 of pressler fc Gunn (1992)| 
but has no spectroscopically measured redshift. There is a published photometric 
redshift ( [I'himm et al 1991 ) of ^ = 0.598±0.11 (galaxy 103) based on nine narrow 
bandpasses spanning 4500A - 8OOOA. This redshift is suspect for several reasons. 
The first is that of the approximately 140 photometric redshifts determined in the 
paper galaxy 103 has the largest quoted uncertainty, two to three times higher 
than than most of the other galaxies. Secondly, the galaxy is classified as type 
Scd based on its spectral energy distribution (SED), however, the luminosity 
distribution in the HST images is inconsistent with a disk system and is an 
excellent fit to a deVaucouleurs r^/^ profile (see below). These two points seem 
to imply that the SED of this galaxy is unusual and since the quoted redshift 
differs by less than 1.5a from the cluster redshift {z = 0.46) we assume the galaxy 
is a cluster member at redshift z = 0.46 for the rest of this paper. 

The position of the lensing galaxy determined using the STSDAS 'metric' 
is (RA,Dec) = (14:11:19.6,52:11:29) (J2000). The best-fit observed parameters 
for the elliptical galaxy are half-light radius Re = 0.61 ± 0.03", ellipticity = 
0.27 ± 0.03, position angle 38.3 ± 1 degrees and total magnitude F702W(AB) = 
20.78 ± 0.05. Using the assumed redshift of z = 0.46 and the [Coleman et al.| 
(1980)1 spectral energy distributions (convolved with a filter response function 



and converted to M^, as in [Lilly et al. 1995]) we can convert to r(gunn). We 
find r(gunn) = F702W(AB)+0.51 for a present-day elliptical spectral energy 
distribution. Thus our F702W(AB) magnitude corresponds to r(gunn)= 21.29 ± 
0.05 which agrees with the previously measured value of pressler fc Gunn (1992) 



of r = 21.20. Their value of Re = 0.6" for the half-light radius also agrees well 
with our value of Re = 0.61" although we find that an exponential disk (their 
profile classification) provides a very poor fit to this galaxy while a deVaucouleurs 
profile fits very well. 

Assuming a redshift oi z = 0.46 leads to Mb(AB) = -19.32 + 51og(/z), Re = 
2.1h~^ kpc, and restframe B central surface brightness /xo,b(AB) = 13.57 mag 
arcsec"^. The adopted k-corrections are for a present-day elliptical galaxy and 
because the observed band corresponds to 4780A in the galaxy's rest-frame (not 
far from B) the effect of a reasonable amount of color evolution will be very small 
(< 0.05 mag). 

An estimate of the expected velocity dispersion can be made by assuming 
that this galaxy lies on the fundamental plane. Surface brightness evolution of 
elliptical galaxies ([Schade et al. 1997]) with redshift will effect this estimate signif- 
icantly. To compare with [Jorgensen et al. (1996)| we first transform the surface 
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brightness to the Gunn r-band. A present-day elliptical has (B-r)(AB)= 1.18 
yielding a central surface brightness /io(r,AB) = 12.39 and r(gunn)=r(AB)+0.21 
yielding /Uo(r) = 12.60. For a deVaucouleurs r^l^ law the surface brightness at i?e 
is /Xe = A^o + 8.33 and the mean surface brightness within R^. is < yUg >= /^e — 1-39 
so that < /ie >r= 19.54 mag/arcsec^. [Jorgensen et al. (1996)| give, for the Gunn 
r-band, log < Jg >= -0.4(< /Xg > -26.4) to yield a value of log < 1^ >= 2.74. 
The fundamental plane we are using is defined at z ~ 0. If we assume a sur- 
face brightness evolution law such as that found in [Schade et al. 19^ then we 
should correct the surface brightness by ~ 0.5 mag to bring it to its present-day 
(non-evolved) value. This yields log < Ie >= 2.54. 

Using our measured value of logi?e = 0.31/i~^ kpc allows us to read an es- 
timate from Figure 1(b) of |Jorgensen et al. (1996]| . If we ignore the evolution 
correction then we obtain 1.24 log a — 0.82 log < Ig >= 0.6 which is equivalent 
to cr = 195 ± 20 km sec~^. On the other hand the evolution corrected surface 
brightness yields a = 150 ± 15 km sec~^. The quoted errors are estimated from 
the scatter in Figure 1(b) of [Jorgensen et al. (1996)| . 

In Table |I] we show the positions of the four lensed background images rel- 
ative to the lens galaxy (0,0). The positions and uncertainties were determined 
by measuring the locations on each of the six exposures and taking the mean 
and standard error. Also given are the F702(AB) magnitudes for each image. 
The uncertainties reflect only photon noise and not systematics due to galaxy 
subtraction errors (which may be significant for image D). The faintest image 
(D) is detected at a signal-to-noise of 12. The galaxy is located in WFC2 and 
the four images appear unresolved. Therefore, the source is likely to be a quasar. 



3. Lens Models 

3.1. Gravitational Lens Theory 

In order to test the lensing hypothesis we have attempted to model the sys- 
tem with simple lens models. To briefly review, for a given source position the 
lensed image positions are given by the extrema of the virtual time delay surface 
( [Blandford fc INarayan 1986| , [Kochanek 19911) : 



At = ^{x-uf + ^{y-vf-(P{x,y), (1) 

where {x, y) and {u, v) are angular coordinates in the image plane and source 
plane, respectively. The two-dimensional effective potential is: 
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DolDosc^ 

where ip is the two-dimensional projected potential and Dls, Dql, and Dqs are 
the angular diameter distances between the lens and source, observer and lens, 
and observer and source, respectively. The extremes of Eqn |1] occur when: 

d(f> , d(p 
u = x—— and v = y — —, [6) 
ox ay 

and the image magnification is given by: 

M-'=fl-^Ul-^)-f^V. (4) 

\ dx"^ J \ dy"^ J \ dxdy J 



3.2. Fitting the Data 

The fitting procedure involved varying lens model parameters and source 
position until the best fit to the measured image positions was obtained. In 
our model fitting we decided to only use the lensed image positions and not the 
relative magnifications due to the high probability of microlensing by stars in the 
lens galaxy ( [Witt et al. 1995| ). The goodnesss of fit was determined from: 

X = 2^ 2 — ' (5) 

errf 

where Hc is the number of constraints, Mj is the predicted value for the observable 
Oi, and ervi is the corresponding uncertainty. In our case = 8 (the four image 
positions). The number of degrees of freedom is DOF = ric — Upar where Upar is 
the number of model parameters (including 2 for the source position). 

The first model we tried was a single lens centered on the galaxy. We assumed 
an edge-on, axisymmetric, oblate, singular isothermal ellipsoid mass distribution. 
The projected mass distribution is ( [Keeton et al. 19971) : 

E(i?, 6) = ~ '''^ (6) 

2GR^4, + l^el, + 1 - 2eM cos2(^^ - 6,) 

where e*/ is the ellipticity (one minus the axis ratio) of the projected mass dis- 
tribution, 9m is the position angle of the projected major axis and Vqq is related 
to the mean line-of-sight velocity dispersion by: 



where a = eA./(2 — €m)/{1 — euY- The corresponding two-dimensional effective 
potential is ( [Kassiola fc Kovner 19911 ): 



cosf 



+ sin(^ - e,) sinh"^ [f^ sin(^ - 9,)] } 



where r = [x^ + y^)^-^ and 9 = ta.n~^{y /x) are the polar coordinates of the 
angular position on the sky with respect to the galaxy center. The parameters 
vgo, Cm, 9e and the source position are varied in the fitting procedure until the 
lens reproduces the observed lensed image positions as closely as possible. The 
best fit parameters along with are shown in Table |^. This model is a very 
poor fit to the data, the value of I DOF is unacceptably high, and the implied 
ellipticity for the galaxy is also higher than one would expect from its optical 
appearance. 

The lensing galaxy is a member of a cluster whose center is located roughly 
perpendicular to the line joining images A and C. The second model we tried 
consisted of two singular isothermal spheres, one centered on the galaxy and one 
centered on the dominant cluster galaxy (offset is 8.66" E and 39.76" N from the 
primarly lensing galaxy) representing the cluster mass distribution. The effective 
potential in this case is: 

'^2 = 4vr^^(4o^ + t;^oO. (9) 

This model produces a better but still very poor fit to the data (x^ / DOF = 
21.19 / 4). 

The galaxy appears elliptical in the WFPC2 image motivating us to try a 
model consisting of a singular isothermal sphere for the cluster and an isothermal 
elliptical potential for the galaxy. 

= + 47r%^r'. (10) 

This model provided a very good fit to the image positions (x^ / DOF = 1.58 
/ 2). Additionally, the inferred values for the line-of-sight velocity dispersions, 
vios = 175 and 135 km s^^ for Zg = I and 3, respectively, are close to the 
predictions from the fundamental plane relationship (see §0). The cluster velocity 
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dispersion is also within the range seen for other similar clusters. However, the 
inferred position angle of the mass distribution major axis differs by about 65 deg 
from the major axis of the luminosity distribution. Attempts to force the mass 
to be oriented identically to the light produced values larger than model 2 
(spherical galaxy) for non-zero ellipticity. Despite this, we conclude that the 
success of the simple double-lens models in fitting the image positions implies a 
high probability that HST 1411.3-1-5211 is a lensed system. 

In fact the lens mass model must be more complex. For example, the individ- 
ual galaxies in the cluster will cause deviations from a smooth mass distribution. 
However, until the redshift of the source is determined (and the redshift of the 
lens is verified) and more constraints arc found (perhaps from VLBI imaging) 
there is little point in adding complexity to the lens model. 

4. Conclusion 

In this letter we report the discovery of a new quadruple lens HST 1411-1-5211 
discovered in deep WFPC2 images of a galaxy cluster a,t z — 0.46. The redshift 
of the primary lensing galaxy is currently unknown but it is quite likely a member 
of the cluster. The images of the source appear unresolved in the WFC implying 
that the source is a quasar. We have successfully applied simple lens models to 
the data. Acceptable fits require both an elliptical galaxy mass distribution and 
an external cluster. Simpler models consisting of only a galaxy or spherical galaxy 
plus cluster do not provide viable fits to the data. The model galaxy velocity 
dispersion is consistent with the inferred value based on the fundamental plane 
relationship, however, the mass distribution appears to be misaligned from the 
light distribution by approximately 65 deg. We conclude that HST 1411-1-5211 
is a very strong candidate for a gravitationally lensed system. Determination of 
the source redshift and confirmation of the lens redshift would allow us to put 
strong constraints on the mass properties of the lensing galaxy. 

Support for this work was provided by NASA through grant # HF-01069.01- 
94A. 
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Fig. 1. — The upper panel shows the combined HST image of the regions around the lensing 
galaxy. The lower panel shows the same regions with a galaxy model subtracted. Each panel 
is 4.3" X 2.0". North is up and east is to the left. 
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Table 1. Image Positions ^ 



Image 


RA 


Dec 


F702(AB) 




(arcsec) 


(arcsec) 


( mag.) 


A 


1.115 ±0.003 


-0.399 ±0.004 


24.96 ±0.03 


B 


0.120 ±0.025 


0.612 ±0.016 


25.00 ±0.03 


C 


-1.125 ±0.007 


0.019 ±0.005 


24.92 ±0.03 


D 


-0.261 ± 0.005 


-0.703 ±0.007 


25.95 ±0.10 



^Galaxy at (0,0) 



Table 2. Parameters and Time Delays for Lens Models ^ 



Model 


u 


V 




vco 


Em 


0. 


A-B 


A-C 


A-D 


/ DOF 




(arcsoc) 


(arcsoc) 


(km s"') 


(km s-i) 




(degrees) 


{h^^ days) 


{h days) 


{h days) 




1 


0.056 


-0.092 


355 (280) 




0.74 


14.9 


171.6 (107.7) 


-4.8 (-3.0) 


-94.9 (-59.6) 


388.40 / 3 


2 


3.567 


16.264 


240 (190) 


1190 (940) 






17.2 (10.8) 


3.8 (2.4) 


9.8 (6.1) 


21.19 / 4 


3 


4.467 


20.441 


245 (195) 


1330 (1055) 


0.45 


103.1 


-4.1 (-2.5) 


-33.4 (-21.0) 


-11.0 (-6.9) 


1.58 / 2 



^Assumes Qq = 1, and Zg = 1. Values in parentheses assume Zg = 3. 



